Longitudinal Modes of Elastic Waves in Isotropic Cylinders 

and Slabs 

By A. N. HOLDEN 

The Rcncral properties of the longitudinal modes in cylinders and sinks nre de- 
veloped with the aid of the close formal analogy between the dispersion equations 
for the two cases. 

1. Introduction 

THE classical exact treatments of the modes of propagation of elastic 
waves in isotropic media having stress-free surfaces but extending 
indefinitely in at least one dimension are those of Rayleigh 1 for semi- 
infinite media bounded by one plane, of Lamb 3 for slabs bounded by two 
parallel planes, and of Pochhammer* for solid cylinders. Rayleigh showed 
that a wave could be propagated without attenuation parallel to the sur- 
face, in which the displacement amplitude of the medium decreased expo* 
nentially with distance from the surface, at a velocity independent of fre- 
quency and somewhat lower than that of either the plane longitudinal or 
plane transverse waves in the infinite medium. Such "Rayleigh surface 
waves" have received application in earthquake theory. 

For slabs or cylinders the treatments lead to a transcendental secular 
equation, establishing a relation (the "geometrical dispersion") between the 
frequency and the phase velocity, which for some time received only asymp- 
totic application in justifying simpler approximate treatments. The past 
decade, however, has seen a revival of interest in the exact results 4 - 4 stimu- 
lated by experimental application of ultrasonic techniques to rods*** and 
slabs/ by the use of rods and the like as acoustic transmission media, and 
perhaps by curiosity as to what qualitative correspondence may exist be- 
tween such waves and the more intensively studied electromagnetic waves 
in wave guides. That this correspondence might not be close could be antici- 
l>ated by observing that an attempt to build up modes by the superposition 
of plane waves in the medium reflected from boundaries would encounter 
an essential difference between the two cases: the elastic medium supports 
plane waves of two types (longitudinal and transverse) with different veloc- 
ities, and reflection from a boundary transforms a wave of either type into 
a mixture of both. 

On grounds both formal and physical it may be expected that solutions 
to the equations of small motion of the medium with a stress-free cylindrical 
boundary can be found with any integral number of diametral nodes of the 
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component ol displacement along the rod, as well as for the "torsional" 
modes in which there is no displacement along the rod whatever. The clas- 
sical results are for no such nodes, the "longitudinal" (or "elongationa!") 
modes, and for one such node, the "flexure!" modes. The secular equation 
for modes with any number of such nodes has been exhibited by Hudson. 4 
For any one of these tyj>es of mode, it may be expected that the secular 
equation will define a many-branched relation* between frequency and 
phase velocity, and that a different number of interior cylindrical nodal 
surfaces for the displacement components might be associated with each 
branch. Apart from the relatively simple torsional modes, the only branches 
whose properties have been intensively studied arc the lowest branch of 
the longitudinal 1 and the lowest of the flexural* modes, because they (and 
the lowest torsional branch) are the only ones extending to zero frequency, 
the others exhibiting "cut-oil" frequencies at which their phase velocities 
become infinite and below which they are rapidly attenuated as they prog- 
ress through the medium. 

Three qualitative results of these studies are of especial interest. In the 
first place, with increasing frequency the phase velocity in the lowest lon- 
gitudinal and flcxural branches approaches the velocity of the Rayleigh 
surface wave, and the disturbance becomes increasingly confined to the 
surface of the cylinder In the second place, the dispersion is not monotonic 
as it ift in the electromagnetic case: the phase velocity exhibits a minimum 
in the lowest longitudinal branch 4 and a maximum in the lowest flexural 
branch* with varying frequency. Finally, in the lowest longitudinal branch 
at least, the cylindrical nodes of the displacement components vary not only 
in radius but even in number with the frequency.* 

The last result suggests that it would be difficult in practice to drive a 
cylindrical rod in that pure mode represented by its lowest longitudinal 
branch over any extended frequency range, since it is difficult to visualize 
a driving mechanism having suitable nodal properties. Longitudinal drivers 
which can be readily constructed may be expected to deliver energy to all 
longitudinal branches, in proportions varying with frequency. How satis- 
factory such a transmission device could be would depend importantly on 
how much the phase velocities at any one frequency differed from branch 
to branch. 

This paper sketches the behavior of the higher longitudinal branches. 
That behavior could, of course, be determined exactly; Hudson* has shown 
how the calculation of the roots of the secular equations can be facilitated, 

* This i* true in particular of the flexural tync of mode, and in his otherwise excellent 
treatment of flexure Hudson's statement to ihc contrary must he disregarded. Recenl 
writings in this field have tended to distingiiMi as "branches" what in allied problem* 
are commonly called "modes". 
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and Huctcr* has used graphical methods* The alternative adopted here is 
a semi-quantitative treatment, assisted by extensive reference to the be- 
havior of longitudinal waves in slabs,* for which the secular equation is 
simpler and closely analogous. The analogy in the case of flexural modes is 
considerably less close and will not be discussed. 

The general consequences of the inquiry are that the higher longitudinal 
branches have phase velocities which are not necessarily monotonic func- 
tions of frequency. With increasing frequency, however, those velocities all 
approach that of the plane transverse wave,** not that of the Rayleigh sur- 
face wave (nor that of the plane longitudinal wave, as some investigators 
had guessed), a fact reflected perhaps in the experimental observation that 
driving a rod transversely usually provides purer transmission than driving it 
longitudinally. t Variation of nodal cylinders in location and number with 
frequency persists in the higher branches. 

2. The Slab 

The slab extends to infinity in the y t z plane and has a thickness 2a in the 
^-direction. The displacements of its parts in the x y y t z directions are u> v t 
tt\ Its material has density p and Lame elastic constants X and >i, so that 
its longitudinal wave velocity is \/(2m + X)/p and its transversa wave veloc- 
ity is VV/p' That /* shou,d be positive is a stability requirement of ener- 
getics; X will also be taken as positive since no material with negative X is 
known* 

The equations of small motion are, in vector form, 

(2* + X) grad div («, v t w) - » curl curl (u t v t to) = p — («, v f w\ 

Solutions representing longitudinal waves propagated in the z direction can 
be of the form 

where U is an odd function, and ir an even function, of x alone, w is the 
frequency in radians ]>cr second, and 7 = <afc where c is the phase velocity. 
Solutions independent of y arc chosen here because they provide the simplest 
analogues to the case of the cylinder. Substitution shows that U - A€ t 

•lam indebted to Dr. W. Shockley for the tufigestion that this behavior might dis- 
play a clo*e enough analogy to that of l he cylinder to provide insight; ihe wort; of Morse 
bear* out the analogy. 

•'The fact is noted by Bancroft. 

1 1'rivate communication from 11. J. McSkimin. 
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W ^ Be*'t is a solution (where A and ft are constants measuring the am- 
plitude), if cither 

(i) *»-*£--./ and yA l -k l B t , 
or (ij) A J - £L - 7 r a nd Mi = ->&• 

M 

When solutions of both types are so superposed as to make (■' odd and H' 

even 

U « Mi sin *i x + iA t sin k t x, (1) 

H' - <4] ? cos k t x - /I, -cos *,*. (2) 



The normal and tangential stresses on planes perpendicular to x 



.ire 



X , = {2M + X)_ + X^+-j, J.-n^ + ^j, 
and the requirement that they vanish at x - ±a leads to the boundary 

conditions 

A x (\y* + (2m + X)if) cr« Aiff + lAipkik* cos **fl = 0, 

2/1, 7* &in *iO+ -4*(r - tj) sin**a = 0, 

the vanishing of whose eliminant with regard to At and A 7 is the secular 
equation. Although in principle that equation establishes a relation be- 
tween 7 and w, it is more conveniently examined when expressed in terms 
of a = ha and (J ■ k^a, which are quadratically related to <j> and y by (i) 
and (ii). In those terms it becomes 

Q&+ (2v+ \)<fy cos a sin 

+ 4(m + AjotffwF - (2/i + \)<r) sin a costf = 0. 

The physically interesting quantities can be expressed in terms of a ai d 
0, with the aid of (i) and (ii). Thus 



"•'-S^w «*-■>• - 



, ftfl* - (2m + X)q' 
f a'Gi + X) 



(4) 
(5) 
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Hence, denoting / ■ 0/a 9 the phase velocity c = <»}/y is given by 

PC = —ft -— - - E t {6) 

HP - (2m + X) 

where E is an "effective stiffness", a function of the elastic constants 
and /. 

Since = is a trivial root of equation (3), it can be divided by 0, and 
the expression on the left then becomes even in both a and and (3) can be 
regarded as an equation in cc and 1 . From (4) and (5) it is evident that, 
if &) and 7 are both to be real, a* and 0* must be real and must obey the ine- 
qualities 

?><*, v& > (2 M + \)<*> (7) 

and thus the root ■■ a can be neglected. The general character of the 
desired roots can consequently be exhibited on a plot of 0* against a*. Evi- 
dently on that plot lines of slope unity are lines of constant frequency 
(equation 4), and lines radiating from the origin arc lines of constant veloc- 
ity (equation 6). As will appear later, however, it is more convenient to use 
a linear rather than a quadratic plot, real a being measured to the right, 
imaginary a to the left, of the vertical axis, and real upward, imaginary 
downward, from the horizontal axis. Here radial lines arc still lines of 

constant velocity, but lines of constant frequency arc no longer simple. 

In Fig* 2 such a plot has been sketched for the first few modes of a mate- 
rial obeying the Cauchy condition X = ;; , the properties shown are restricted 
to those derived in the following paragraphs, and are lettered in Fig. 1 to 
correspond with those paragraphs* 

(a) By virtue of (7), the significant portions of the roots lie above and 
to the left of the lines 1 - a\ tf* - (2m + X)rf. Setting «F - (2ji + X)a* 
in (3) reveals the cut-offs at sin = and at cos a — 0: in other words at 

£• = nV, a 1 - ~r-i «***> and also at fl 3 - 2 ^ + A 
2p + X M 




■(■ 



IV 
+ ~ J t 3 , where n is any integer. 



(b) Setting a - in (3), it can be seen that the roots intersect the line 
a* — at the points sin — 0. By calculating the derivative of 1 with re- 
spect to a 3 , those points (at which a changes from pure real to pure imaginary) 

can be shown not to be multiple points, and the branches to have — = 
and -jj-j, = — v t y independent of branch number and negative for 
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Fig. 1— Lines and intersections, discussed in the correspondingly lettered paraicraph* 
ol the text, which determine the properties of the first five branches of the lonRiiudn al 
modes for a materia) obeying the Cauchv condition. Two coincident p«in* of (Mints are 
marked (2). 
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Fir. 2 — A rough sketch of the branches determined by the properties illustrated in 
Fig. 1. 
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all materials* At those points the phase velocity is that of a plane longi- 
tudinal wave. 

(c) Setting sin ■ reduces (3) to afl(w9* - (2m + XV) sin a - 0, 
and hence in the region of positive 0* and negative a 7 the roots do not inter- 
sect the horizontal lines 0* = n V (ft ^ 0). As can be seen from (6) this 
confinement implies that the velocity is asymptotic to Vm/p* that °' a 
plane transverse wave, with increasing frequency, Notice that, in the re- 
gion of ]x>sitive ff and negative a* t takes the values cos - only where 
\fi* + (2^ + XV = and that the roots have zero dope there. At those 
|x>ints the waves have a phase velocity V2 times that of a plane transverse 
wave, 

(d) In the region of positive 0* and positive a* the roots exhibit a some- 
what more complicated behavior, but confining lines can again be found: 
the diagonal lines S - (« + i> - a- It is the nature of such critical lines 
as these which can be better exhibited on the linear than on the quadratic 
plot. Alternatively those lines can be written cos a cos — sin a sin ^ 
(and thus will be shown lo have analogues in the case of the cylinder), and 
substitution of this expression into (3) shows that if the roots intersect 
these lines they must do so for values of a and satisfying the relation 

4(m + X)a0(w3* - (2m + XV) = -(X0* + (2m + XV) 1 cot* a. 

Hut the inequalities (7) make such values impossible, 

(e) This suggests that in that region the roots may oscillate in a some- 
what irregular manner about the diagonal lines = «x — cr. Indeed it is 
immediately evident that they pass through the points cos = cos a = 
and sin — sin a = 0. 

(f) Expressing those lines as sin a cos + cos a sin 0=0, and sub- 
stituting into (3), shows that additional intersections may be afforded by 
any roots of the quarlic equation 

OyP + (2m + XV)' - 4(m + \)a»(i& - (2/i + XV) - o 

which obey the inequalities (7). Discarding the root a + = 0, and dividing 
by a 3 , yields the cubic equation 

XV - (2 M + X)V + (2m + X)(2m + 3\)t + (2m + X)* = 0, (8) 

whose roots are the negatives of the roots of the cubic equation for the 
Raylcigh surface wave velocity- It is well known that the Raylcigh cubic 
always yields one and only one significant positive root, and hence equation 
(8) can afford at most two additional significant intersections of any root 
of the secular equation with the line about which it oscillates. Although it 
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is not feasible to exhibit the roots of the Rayleigh cubic explicitly for arbi- 
trary u and \> it is of some interest to exhibit its discriminant 

D - 1 \\2n + X) 3 (m + X)'(11A 3 + 4AV - 9Am* - 10m 1 ), 

and to note that for real positive values €>f X and ji it changes sign only once, 
at approximately X/m ■■ 10/9. Hence for X//i > 10/9 two roots of the Ray- 
leigh cubic arc complex, while for X/ji < 10/9 two roots are real and nega- 
tive. For a material olwying the Cauchy condition X/m — 1< the roots of 
the Rayleigh cubic arc -3, -3 ± 2V^; thus I - 3, 3 + 2\/3, both of 
which ol>ey the inequalities (7), are relevant to intersections of each branch 
of the roots of the secular equation with the line about which it oscillates, 
(g) The results of (e) and (f) suggest the value of a similar investigation 
in the region of imaginary a. Here (denoting a = iA and L m ft/A where 
A is taken positive and real) intersections occur between the branches and 
the lines sinh A cos - cosh A sin & - when (XL* - (2ji + X)) 1 = 
4fr* + X)/,Gi£* + (2/i + X)). Clearly this quartic in L has two an d only 
two positive real roots, one greater and one less than y/(2$i + X)/X, In the 
case X = ji, those roots arc approximately 9 and 1/3. This information, 
taken with that of (c), establishes that the branches arc confined in the 
region of imaginary a to bands determined by nv < < (« + j)w, having 
one tangency to the lines = (« + $>r; and that at values of A greater 
than correspond to the smaller root of /, the branches lie in the hands ht < 

& < (« + frr. 

It is convenient to obtain assurance that in general the branches do not 

intersect at any point by noting that the confining lines of paragraphs (c) 
anil (d) define bands within each of which in general one and only one cut- 
off point falls. Pivoting a ruler about the origin of Fig. 1, and recalling the 
cut-off conditions, avails* Degenerate cases arise when the elastic constants 
satisfy a condition 2ji + X - ntyi where ft is an integer; in those cases some 
cut-off points coincide in pairs on some confining lines. Calculation of de- 
rivatives at those jioints shows that the cases are not otherwise exceptional: 
the pair of roots forms a continuous curve which is tangent to the cut-off 
line at the double cut-off point. 
From (6) it follows that the phase velocity will have a maximum or a 

minimum with frequency if-—; ■= -=- That condition requires 

d\a') or 

, (X fl* + (2u + X)aV(xV + (2g + X)(2m + 3 X)*") 

40* + X)(2^ + X)VGP - rfJW 1 - (2*i + XV) " 



OW 
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In the region of positive 0* and a* and of the inequalities (7), this is impos- 
sible, but when a* is negative the condition may be satisfied. If it is satis- 
fied in the higher branches, however, it must be satisfied an even number 
of times m any branch, so that the branch exhibits as many maxima as 
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Fig. 3— The beginnings of the dispersion turves inferred from I*ir> I and 2. Tie soid 
ahaded line* are the curves about which the branches oscillate, interceding them it the 
triangles and lyin£ on the shaded aide of them elsewhere, and the dashed extensions are 
the branches themselves. For increasing r these branches aU become asymptotic to the 
base line. The dashed Lines at the top are the true cutoff frequencies; the solid cut-off 
lines arc the asymptotes of the shaded curves. The bc^innm* of the lowest branch is 
shown at the lower left; it becomes asymptotic to a line below this pk>t t 



i-\Z*-vi- 



ifter passing through a shallow minimum. 



minima, for clearly the phase velocity is a decreasing function of frequency 
near the cut-off, and the velocity can also be shown to approach its asymp- 
totic value at high frequencies from above in the higher branches. 

In finally displaying the dispersion curves (Fig, 3) it is convenient to use 
as reduced variables r t the number of plane transverse wave lengths in 
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one slab thickness (which is proportional to the frequency), and **, the 
ratio of the velocity to that of a plane transverse wave. Evidently 
* == A*Y .. l 2m + X (J ~*V ( C Y « (2m + X)(/ j - 1) 

W ** 7T\ {tr '" h W & - (2m + x) ■ 

where c* s Vm/p- 

For completeness, the lowest branch will be briefly sketched: that which 

originates at a = — 0. A calculation of ^ at that point yields only one 

non-trivial root, -(2 ** + X)(2m + 3X)/X*, and thus the phase velocity at 
low frequencies is found to correspond, as would be expected, with that 
given by the stiffness (a semi- Young's modulus, so to speak) of a material 
displacement-free in the ^-direction but not in the y -direction, E = 
4jiG* + X)/(2 m + X)* Since lines radiating from the origin of the (0*1 a 1 ) 
plot arc lines of constant velocity, the dispersion curve for this branch 
starts with zero slope. The root curves over, intersecting the line 

X0* + (2m + X)a* - at & = ? t and intersects die line & - again at 

-4* (or ■ iA) where (2/i + X) A cosh A = 4(m + X) sinh A. For large 
negative «* and fl 1 , equation (3) approaches 

XV* + 4^0* + X)P + 2X(2^ + X)J* - 4(2p + X)G* + X)J + (2 M + X)* = 0, 

which after discarding the trivial root I = 1 leaves the Raylcigh cubic. In 
the case of this one branch, the phase velocity approaches its asymptotic 
value at high frequencies from below, and hence the dispersion curve must 
have an odd number of maxima and minima, and in particular at least one 
minimum, as was discovered by numerical calculation for the corresponding 
branch in the case of the cylinder by Bancroft, 4 

The complicated behavior of the displacements in the higher branches is 
sufficiently illustrated by a brief consideration of their nodes: the values of 
x at which the displacement is entirely along or entirely across the slab. 
From (1) and the boundary conditions, the ^-dependence U of the displace* 
ment component perpendicular to the slab will be given by 

K % U = (X0* + (2 M + X)a a ) sin sin — 

a 

(la) 

+ 2(jt& - (2^ + X)a*)sin a sin - 
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or by 



ax 



K%U = 2(n + X)a0 cos sin — 



- (XJJ* + (2*» + X)o*) cos a sin ^ 



(lb) 



where AT, - (X^ + (2m + XV)^ sin ft A", - 2(/x + XjaflAT cos ft at all 
points (o, (?) satisfying (3). Similarly from (2) and the boundary conditions, 
the x-dependence W of the displacement component along the slab will be 
given by 

K, W = <X0* + (2 M + \W) sin /J cos ^ - 2(m + X)<tf sin a cos ^ (2a) 

or by 

K t W = 2(jt? - (2m + X)a') cos cos — 



(X0 1 + (2 M + X)a*) cos o cos ^, 






(2b) 



A", = ia>o(M + M ^ + ^t^ A-sinftg t - 2iara0i + X)tf cosft 

Examine now, for example, a material for which X = ji, in the branch 
whose cut-off is at fi « 2x t a - 2v/y/3. It can be verified at once that the 
nodes of the two components at some of the values of {a, &) discussed ear- 
lier arc described by the following table f in which / = cos — J : 



a 





2-T/V3 


2r 


A' 


2r 


wfl 


Sw/2 





3t 

7,/2 

3r 



Value* of */* for aode* of tf 



0, ± two values given by 

2/V3 sin 2wx/aV$ - 

cos 2wfV3 sin 2«/a 
f ±1, ± one value given by 

141 + 8 — 
0, ± two values fliveii by 22/* + 

11/- 2 - 
t ±J,±|, ±\ 

0, ±i, ±| ( ±1 



Value* of i/j (or oodrt of ff 



±i.±i 



± two values given by 

I4/ 1 - 2/ - 7 - 
±l» ± one value given by 

5J 1 - 15/ - 11 - 
no nodes 



It is to be noted in general that the nodal variations become less extreme 
at high frequencies, since for all branches except the lowest U and W tend 
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lo become proportional to sin - and cos — respectively, and £ approaches 

the value nr where u is the branch number in order of increasing cut-off 
frequency, with « = ascribed to the lowest branch. Thus the feeling, de- 
rived from more familiar cases of wave motion, that the order in which 
the branches arrange themselves should be correlated with the number of 
nodes they display retains an asymptotic validity here, in respect of each 
displacement component. 
Xodes of absolute displacement will occur only at special frequencies. 

With the notation a' = t ff ■ — , the conditions for their occurrence can 

a a 

be written 

(tfj* - (2? + X)a") sin & cos <J + fci + X)aV cos & sin a' => 0, 
sin 2a f sin Iff 0' & # - A , . 



taken together with (3). 



J. Tut; Cvlinpik 



Procedures analogous to those of the preceding section, and presented by 
1-ovc 1 , lead to Pochhammcr's secular equation, which in the present 
notation* is 

(XJF + (2ji + \)c??Ma)JiW 

+ 4(y + X)o0(«3» - (2*. + X)o*)/,(a)y.C8) (9) 

+ 20* + X)(2m + X,)a(«» - &H«)Jxifi) - 0, 

where (4), (5) and (6) still hold, with o signifying the radus of the rod The 
analogy between ($) and ihe first two terms of (9) is striking- Again the 
roots = and = a can be neglected, and the equation when divided 
by /3 becomes even in a and &* and a plot of ff 1 against cr becomes appro- 
priate, with the restrictions (7) as to regions of significance. The following 
paragraphs are lettered to correspond with their analogues of the preceding 
section. 

(a) Setting tfp = (2/i + \W in (9) reveals the cut-offs at Ji(fi) = 
and at (2 M + \)*J (a) - Wito- 

(b) Setting it - in (9), it can l>e seen that the roots intersect the line 

• In comparing this treatment with that of Hud*on\ interpret his symixris 
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a 1 - at the points /i(0) = 0, traversing them with j^l - — ^ ', 

or half that of their analogues. Again it is only at those points that the 
phase velocity has the value VT2jt + ^)/p* 

(c) Setting /,(0) = reduces (9) to apbF - (2p + \W)M<x) - t and 

hence the roots arc confined between the horizontal lines J\(fi) ^ in the 
region of positive 1 and negative a*, and the velocity is asymptotic to 
VV/p w ' t ' 1 increasing frequency. They meet the line XJ3 3 + (2j* + X)o* ■ 
at the points &Jolfi) — Jiffi) ■ 0, or in other words at the maxima and 
minima of /if/?),* where again the phase velocity is V%/p- 

(d) In the case of the cylinder confining lines** are 



[ 



*\Nr + (2/* + X) a 1 ) 



J[(fi) - AMAifi) = 0. 



since substitution shows that intersection of (9) with these lines would re- 
quire 

(Xf + (2*. + X)a*)» " 

L a(X^ + (2w + X)a') J 

which cannot be satisfied by permitted values of cr 1 and &. 

(e) This suggests that in that region the roots may oscillate about the 
lines 



JM/M + JM [/.W + ^f (i> ;\ M «<■> 



a 



In fact, in view of the equivalence j\(x) = /o(x) — - Jx(x), those lines 

ran be seen to have points in common with the roots of the secular equation 
at /i(a) = 0, Jilfi) = 0, and at 

fM-o. JIM + ^^^JM - a 

(0 Substituting the expression for the lines of (e) into (9) shows that 
again additional intersections may be afforded by suitable roots of the 
cubic equation (8). 

'The analogy to the corresponding intersections (or Lhe slab at ihe maxima and 
minima of sin is noted by Lamb. ref. 2, p. 122* footnote. 

** There arc infinitely many such families of lines but none carries the analogy with the 
slab to the point o f being independent of the elastic constants. The families used in (d) 
and (e) serve the present purpose as simply as any. 
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(g) In the region a — iA t the analogous lines arc given by 

(MUVSU + /.<*> [liiu) - ^tf^ wo] - o, 

with which intersections occur for the same values of &M as in the slab. 

These results permit visualization of a counterpart to Fig, 1 for the cylin- 
drical case. In it the critical lines radiating from the origin are the same. 
The horizontal lincs p instead of being evenly spaced by x/2, arc spaced as 
the zeros, maxima, and minima of /i0). The vertical lines, again no longer 
evenly spaced, are replaced alternately by straight vertical lines J\(a) = 

2(m + X)a 
and by the curved "vertical" lines /*(«) = rip"; — 7; — . \;',/i(a) which 

Ap *f* (Zji x AJor 

lie between their straight companions and approach /«(a) ** as & be- 
comes large. Finally the confining lines, and the lines about which the 
branches oscillate, become the curves defined in (d) and (e), which can be 
seen to follow a course not dissimilar to the diagonal course of their prede- 
cessors, passing through the intersections of the new horizontals and 
verticals. 

Again the branches do not intersect, except for )>air-wise coincidence of 
cut-offs on one or another of the lines (d) when the elastic constants obey 
special relations. Thus the dispersion curves to which Hudson* assigns 

certain of Shear and Fockc's data cannot be taken (and indeed Hudson 
docs not suggest that they must be taken) as corresjxjnding to higher 
branches of the longitudinal modes, since the former curves intersect one 
another, and the latter cannot unless anisotropy modifies their behavior 
qualitatively. The assignments could represent modes other than longitu- 
dinal. The more recent results of Hueter* show essentially the behavior of 
Fig. 3. 

In view of the closeness of the analogy thus revealed, it may be taken as 
probable that qualitative correspondence will obtain quite generally be- 
tween the longitudinal modes of the slab and those of the cylinder 
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Frequency Dependence of Elastic Constants and Losses 

in Nickel 

By K M. BOZORTH. W. P. MASON and H. J. McSKIMIN 

The clastic constants o( nickel crystals and their variation with magnetic field 
(A£ t i : vi r . have been measured by a 10- megacycle ultrasonic pulsing method. 
Tlie con st an is of three crystals agree well with one another when the crystals are 
magnetically aaturaled, but vary with domain distribution when demagnetized. 
The maximum d£ efteel observed is much less (3%1 than has been observed al 
lower frequencies (20%). By measuring the &E effect and the decrement of poly* 
crystalline rods at low frequencies, it is shown that the small effect observed al 
10 megacycles is due to a relaxation in the domain wall motion due to micro- 
eddy-current damping. 

From the initial alope of the decrement- frequency curve, and also from the 
frequency of maximum decrement, the sice of the average domain is found to be 
altnut 0.O4 mm. Actual domains in single nickel crystals have been ohserved 
optically by Williams, who finds domain widths of 0*02 to 0*2 mm. 

THE three elastic constants of nickel have been delermined in several 
single crystals by measuring the velocity of pulses of elastic waves of 
frequency 10 mc/s and duration 0.(X)1 sec. The method has been described 
by McSkimin 1 and the preliminary results on nickel have already I wen re- 
ported briefly,' 

It is well known that Young's modulus, E> increases with magnetization, 
and changes in E (the "&E effect") by 15 to Mi |>cr cent have been observed 
at mom temperature ami changes by greater amounts al higher tempera- 
tures. 3 It was surprising to find then, in our own experiments al 10 mc, that 
the greatest change was only about 3 per cent. It then occurred to us that, 
at such a high frequency, relaxation of the domain wall motion by micro- 
eddy-current flumping might be expected. This led to the investigation of 
the frequency dependence of £E and of the logarithmic decrement, i, in 
polycrystalline nickel, and the results obtained support the theory and give 
information about domain size, as described below. Calculations* based on 
the equations of domain wall motion give results which agree with the ex- 
periments, 

A number of experiments 3 have already established the existence of micro- 
eddy-current losses in magnetic materials subjected to elastic vibrations. 
These losses have their origin in the local stress-induced changes of mag- 
netization of the domains of which magnetic materials are composed. The 
change in magnetization of one domain -will give rise to eddy -currents around 
it and in it, and the consequent loss in energy depends on the frequency/ 
and the resistivity R t and on the size and shape of the region in which the 
change in magnetization occurs. These losses are in addition to the macro* 
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eddy-current losses, due to the relatively uniform changes in magnetization 
of a magnetized specimen that occur during a change in stress. 

Calculations of the logarithmic decrement, 4* , attributable to micro- 
eddy-currcnts, have been made by Becker and Doling 9 and by one of the 
writers.* According to these calculations when the material Ls composed of 
plale-likc domains of thickness/, in which magnetization changes by bound- 
ary displacement, the decrement for nickel, which has its directions of easy 
magnetization parallel to [111] directions, is given by the relation 

5/: ' Uh - 



r» + 3cm. 1 + Pffi 



(1) 



where /o, the relaxation frequency for domain wall motion is/# = ~ — ~' 

/j is the saturation magnetization, E* is the saturated value of Young's 
modulus, mj is the initial permeability, R the electrical resistivity, Am the 
saturation magnetostriction along the (1 11] direction, and c l{ , c t % and c 4 4 
the three elastic constants of nickel which arc evaluated in this paper. For 
low frequencies the initial slope of the -decrement vs frequency curve is 



& . 24£,n3/*X?i, 



j 



5CW (2) 



/ 5«*/J .c„ - c e + 3<r« 

As the frequency is increased the decrement rises to a maximum and then 
declines asymptotically to zero. Both the initial slops of the 5 vs/ curve and 
the frequency at which the maximum occurs are measures of the domain size. 
The initial slope has already been used to evaluate the size of the domains 
infiS Permalloy. 6 It is shown in the present work that the maximum occurs 
in polycrystalline nickel at a frequency consistent with the dimensions of 
domains observed by Williams and Walker* in single crystals of nickel. 

Elastic Constants and Damping is Sikcle Crystals 

The nickel crystals used here were grown by slow cooling of the melt in a 
molybdenum wound resistance furnace, by a method previously described. 7 
They were cut with major surfaces parallel to (110) planes and were placed 
between two fused quartz rods as shown in Fig, 1. Measurements of the 
elastic constants were made as described in detail by McSkimin, 1 by meas- 
uring the velocity of propagation of 10 mc pulses. In order to obtain a num- 
ber of roil actions in the crystal, films of polystyrene approximately J wave- 
length thick are placed between the rods and the nickel crystal. This has the 
effect of lowering the impedances next to the nickel to small values and 
hence nearly |>erfect reflections at the two surfaces arc obtained. The fre- 
quency is varied until successive reflections occur in phase, and the velocity 
is then calculated from the frequency and the dimensions of the crystal. 
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Fig. 1— Experimental arrangement for determining the clastic constants and i£ 
effect in single nickel crystals, 
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The velocities for one demagnetized crystal 1 were found to have the values 
shown by Table I These values of velocity, and a density of 8-90 for the 
single crystal, give values for the demagnetized clastic constants of 

Cu = 2.50, 4s = 1*60- *i = 1-185 (3) 

all in 10" dynes/cm 1 . 

To obtain the AE effect, the whole unit was placed between the jaws of a 
large electromagnet. Since the crystal was about 2.5 centimeters in diam- 
eter but only 0.472 cm thick, saturation could be obtained more easily 
along the long directions of the crystal. Figure 2 shows the changes in 
velocity of propagation along the [1 10) direction, caused by magnetization 
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along |001), for the shear mode with particle motion along [110]. Fields of 
tbotlt 10,000 were attainable but a maximum field of about 6000 was usu- 
ally used. The velocity increases by 2.6 per cent at the saturated value. 
On decreasing the field to zero the velocity drops below the original value 
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Fig. 2— Change in velocity in percent from demagnetized value as a function of the 
magnetizing field for a shear wave in a (110) sscbon when ihe particle velocity t is 
along the llTOl direction and the field U along the 1091) direction* 



for the demagnetized state, but it has practically the initial value when the 
crystal is again demagnetized. The lower value of velocity for the return 
curve indicates that the free energy is lower for some arrangement of the 
elementary domains other than the demagnetized state. 

Figure 3 shows the attenuation in decibels per trip as a function of mag- 
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nctization. The loss drops from about 7 ilb to 1 rib as the crystal becomes 
magnetised. The low value is the remanent loss caused by the energy 
lost to the terminations, so that one can say that the losses due to micro- 
eddy-current and micro- hysteresis are 6 db per trip or 12-7 dbper centimeter 
for this mode of motion. The (> of the crystal can be shown to be equal to 
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Fig. 3— Loss jier Irip (0*472 cm) as a function of magnetizing field for shear wave 
of Fig. 2. 



the phase shift in radians divided by twice the attenuation in nepers per 
cm, or 



2wf/v 



2 (db per cm)/8.68 



(4) 



and our results give Q — 94, corresponding to a decrement of w/Q m 0.033. 

Figure 4 shows a measurement of the same mode when the field is applied 

along the [110J direction. The velocity approaches a slightly different 
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limit on account of the "morphic* 1 effect discussed in another piper 1 If 
we average the two values the effective elastic constant for saturation be- 
comes ¥ 

cU - rfi = 0.954 X I0 l * dynes/cm 1 . (5) 

Measurements for the field along the thickness did not produce saturation 
an:| are not shown, 
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Figures 5 and 6 show similar measurements for the other shear mode 
(Shear 2 of Table I) for two directions of the magnetic field. Averaging 
the two limiting values, the constant Cu at saturation becomes 

cU = 1.22 X 10" dynes/cm* (6) 

The {) and decrement for this case Iwcomc approximately 110 and 0.02S. 
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Figures 7 and 8 show measurements for the longitudinal mode. Variations 
of about 0.6 per cent in the velocity are obtained, and the saturated elastic 
constants, Q and decrement are 

clx + c{t + 2fili - 6,55 X 10 l *. Q - 390, 5 = 0,008 (7) 
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Combining the elastic constants, the saturated elastic constants are evalu- 
ated: 

dx= 2.53, cit - 1-58, cU=* 1-22, (8) 

all in 10 u dynes/cm*. 
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It is obvious from the measurements of Figs. 4 to 8 that the changes in 
the elastic constants with magnetization are much smaller at the high fre- 
quencies (10 mc) than they are at the lower frequencies of 10 to 50 kilo- 
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cycles where changes of 15 to 30 per cent have been observed in polycrystal- 
line material. 3 A rough comparison of the low-frequency values with the 10 
megacycle values can be obtained if we convert the observed changes in the 
c's to the equivalent change in £. This can be done if we use the method 
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developed by one of the writers 9 for obtaining the clastic constants of a 
polycrystalline rod from the cubic clastic constants. In this case the Lame 
elastic constants are given by the formulas 



3 2 4 

+ 2m = j cii + ^ *ii + 5 c M , 

3 ( Cn — c u 
~S C « + — I" 



(9) 
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Since in terms of the Lam£ clastic constants, the value of Youngs modulus 



is 






V 



(in) 



one finds that the difference lietween the saturated and demagnetized value 
of Young's modulus divided by the demagnetized value is 
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AE E^_& 2,381 - 2,312 
£a £» 2.312 



= 0.03 = 3% 



(11) 



which is much smaller than that given by low frequency measurements. 
To check our results, and be sure that the crystals were free from imper- 
fections and strains, two other crystals were prepared and carefully an- 
nealed at 1 100°C The values found for the changes in elastic constants were 
considerably teas for these crystals. The Q f * of the crystal were also higher. 
Table II shows the measured values and the equivalent AE/E values. The 
table shows also the measurements for the demagnetized crystal of two 
Japanese workers 10 1l and the equivalent AE/E assuming that the saturated 
elastic constants are the same as those found for the other three crystals. 
Since these vary by only ±0.5 per cent among themselves, this appears to 
be a good approximation. 

Table II 
Elastic Constants (in 10" dykes/cu*) and a£ Errecr in Single Crystals or Nickel 





UajnctkaUy 5atur*Ud 


1 ■ : i L . : . i li Jed 


». p\ 


CrytUl 










aB/a 


(u 


<U 


In 


■u 


cm 


c« 




1 


2. 53 


1.53 


1.22 


2.50 


1.60 


1.185 


003 


2 


2.524 


1.538 


1.23 


2.52 


1.54 


1.229 0.0017 


3 


2.523 


1.566 


1 -M 


2.517 


1.574 


1.226 0.0046 


Yamamoto" 








2.44 


1.58 


1.02 


0.16 


Honda and 








2.52 


1.51 


1.04 


0.11 


Shimkawi" 

















The lower values of AE/E for the second and third crystals are probably 
due to larger domain sizes, caused by the longer anneal. 

Damping and A£-Effect in Polycrystalline Rods 

To test the theory of micro-eddy-current shielding (see Introduction), 
the velocity and attenuation of elastic vibrations in well-annealed poly- 
crystalline nickel rods were measured over the frequency range of 5 kilo- 
cycles to 150 kilocycles. In the method of measurement/ 1 shown by Fig. 
9, two matched piezoelectric crystals of resonance frequency corresponding 
to integral half wavelcnglhs along the rod, arc attached to the ends of the 
rod. Phase-amplitude balance was obtained by critical adjustment of fre- 
quency and output of the calibrated attenuator. The corresponding level 
was then compared with that obtained when the two crystals were cemented 
directly together. With little error, the velocity of propagation is given by 

_2JU 



»- 1,2,3 



(12) 
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The attenuation A (and hence the Q of the rod) was obtained by solving 
the equation 



s;nh -4/# ■■ 






(13) 



in which r is the ratio of output with crystals together to output with 
specimen attached, U is the length of rod, M c the mass of either crystal, 
Mr the mass of the rod, and Q c the Q of the crystal as determined by reso- 
nance response method. 

For this equation to apply accurately, the terminating impedance pre- 
sented to the rod by the crystals at resonance must be small compared to 
the characteristic impedance of the rod, and the Q of the rod should be >10. 
This method may be used even when the total loss in the rod is so high that 
well defined resonances no longer exist. At the lower frequencies, however, 
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TO 
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Fig. 9 — Experimental arrangement for measuring Che &E effect and associated loss in 
a poly crystalline rod at low frequencies. 



a useful check may be made by the resonance response method of determin- 
ing Q which involves determining the frequency separation A/ for two fre- 
quencies 3 db from the maximum response frequency, and using the formula 



/. 



Ma* 



*/ 



(H) 



Correction for the mass and dissipation of the piezoelectric crystals must of 
course be made. Both methods have been found to agree within about 10% 

the probable error to be expected. 

The Appendix lists formulae to be used when the resonance frequency of 
the crystal driver differs from the frequency at which phase balance is ob- 
tained. This condition of necessity occurred when the rods were subjected 
to a magnetic field, which caused an increase in the velocity of propagation. 

Figure 10 shows a typical measurement of change in frequency and 
change in decrement with magnetizing lield excited in a solenoid surround- 
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ing the nickel rod. Saturation is not quite obtained so that the d£ effect 
measured is slightly lower than the true value, but for relative frequency 
comparisons this is not important. 

The firsl rod measured was 0.320 cm in diameter and 10.16 cm long and 
was annealed at tl00*(\ Five frequencies ranging from 22.5 kilocycles were 
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used and the ratio of the change in Young's modulus to the value of Young's 
modulus for the demagnetized rod is shown plotted in Fig, II* This figure 
shows also the decrement & = */Q. I"t is obvious that the decrement even- 
tually decreases as the frequency rises, and this is contrary to the simple 
theory of the micro-eddy-current effect,* which indicates that the decrement 
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should increase linearly with the frequency* The indicated maximum for this 
rod is below 120 kilocycles. 

In order to obtain the first part of the decrement vs frequency curve, a 
rod of 46.05 cm length and 0.637 cm diameter was next used. This rod was 
annealed at 10S0°C and presumably has a smaller average domain size than 
the first one, so that the important variations occur in a more favorable 
frequency range. 

The changes in elastic constant and the decrement for this rod arc shown 
by Fig. 12 for frequencies from 5 kilocycles to 96 kilocycles. At the lower 
frequencies the decrement increases in proportion to the frequency in 
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Fir. 12— Fractional chanjjc in Young's modulus, and the decrement, plotted as a 
function of frequency for rod No. 2. 



agreement with the simple theory. By extending this curve down to zero 
frequency it is seen that a micro-hysteresis effect (which is independent of 
the frequency) gives an initial decrement of about 0.010, The decrement 
rises to an indicated maximum at somewhat more than UK) kilocycles and 
the change in elastic constant with saturation decreases with frequency. 
The data on these two rods taken together indicate that there is a fre- 
quency of maximum decrement and for frequencies above and below this 
the decrement is smaller. The AE change in the elastic constant decreases 
as lite frequency increases and for very high frequencies the AE effect be- 
comes very small. As shown by the discussion in the next section, the frc- 
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quency of the maximum value of & and the initial slope of the decrement 
frequency curve arc connected with definite domain sizes which can be 
calculated approximately and compared with magnetic domain powder 
patterns. 

Discussion 

Our determinations of the elastic constants may be discussed in relation to 
the values obtained by others. The results reported by Honda and Shira- 
kawa 10 and Yamamoto 11 were unknown to us and unavailable at the time of 
our preliminary communication. The data of the Japanese, converted from 
j-constants to ^-constants by the relations: 

*ii + Jii 
Sn + $nSi% — 2su 

= ~ Jt1 (15) 

Cm - !/*« 

are included with our data in Table II* 

As our experiments show, the 10 mc pulses that we used arc so rapid that 
micro-eddy- currents largely prevent the stress- induced changes in mag- 
netization from penetrating the domains. Therefore the constants deter- 
mined by this method are those for material almost saturated. The values 
at saturation are independent of the initial domain distribution, and of the 
ease with which the magnetization in the separate domains can be changed 
by stress, consequently they are the more fundamental elastic constants of 
the material. The variety of values for unmagnetized nickel is made evident 
from the scatter in the ratios of AE/B lhat have been reported.* The varia- 
tion in the values of the c-constants recently published is thus not surprising. 
The values at saturation of the three crystals examined by us are in sub- 
stantial agreement, as shown in Table II, They cannot be compared directly 
with the results of the Japanese workers because the latter reported data 
for unmagnetized crystals only and then E is sensitive to heat treatment 
and domain configuration. 

As mentioned in the introduction, the damping of clastic vibrations by 
micro-eddy-currents is proportional to the frequency at low frequencies 
(Eq. 1) and it rises with frequency to a maximum and then declines toward 
zero* The frequency at which the maximum occurs has been calculated* 
by using the equation of domain wall motion and evaluating the constants 
from the initial permeability and the power loss caused by domain wall 
motion. The maximum value of & comes at the same value as that calculated 
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for eddy current tosses in sheets having the same thickness as the domain, 
namely 

Pmo/-/* - 0,13 (16) 

/ being the thickness of the sheet, mo the initial permeability, and R the 
resistivity of the material (all in c.g.s. units). 

As noted below, the domain sizes calculated from the initial slope of the 
5 vs/ curve of Fig, 12, and from the frequency at which the maximum decre- 
ment occurs, are respectively 0.035 mm and 0-045 mm (for plates). These 
values agree quite well. The decrement curve is broader than would be cal- 
culated from equation (1) for a single domain size. This agrees with the 
optical measurements of domain size by Williams,* which arc shown by 
Fig* 13. This indicates domain sizes from 0,01 to 0.2 mm. 

The maximum value for the decrement calculated from equation (1), 
using the measured values, is 0.35 compared to the observed value of 0.11. 
Part of this is due to the broadening of the peak caused by a distribution 
of domain sizes, but part may also be Hue to the deviation of the actual 
domain shape from a sheet which has been assume*! in making the cal- 
culations. 

The calculations of domain size arc made in more detail as follows: 

According to Dtiring 13 the change in Young's modulus for nickel contain- 
ing only small internal strains is related to the initial permeability, /i ?1 
as follows: 



Afi XmU - 1>£. 5c« 



j 



£. Swl\ Lfe " <" + 3c. ' (18) 

provided the averaging over all crystallites is carried out with constant 
strain. (If constant stress is assumed, the fraction in brackets, equal to 
1.76. is omitted.) For nickel Xiu is 25 X Ifr*, h is 484, and the c's arc the 
elastic constants given in Table I. This equation holds for low frequencies 
at which the shielding in single domains is negligibly small. When the re- 
laxation effect of domain wall motion is considered 4 equation (18) has to be 
multiplied by the factor 

The data of Fig. 12 give the values: 

£o= 1.83X10" £.- 2.22 X 1 0" dynes/cm*, ^ = 0.21 (20) 

for low frequencies. Using these in the above equation, the calculated value 
of ai„ is 320. A direct measurement* of m has been made for this rod and 
found to be 340, in good agreement with that deduced from the AK effect. 

* Measurements were made independently by Miss M. Goertz and Mr. P. P. Ciofli 
in order to check this unusually high value* 
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Since the permeability is much higher than can be accounted for by 
domain rotation it is obvious that domain boundary motion is occurring. 
Hence in determining the domain sizes from the slope of the decrement vs 




Fig. 13 — Photograph of domains in a single nickel crystal (after William*). Field of 
view, 0.5 mm. 



frequency curve, equations (1) and (2) for domain boundary motion arc 
appropriate. 

When the data of Fig. 12 are extrapolated to zero frequency it appears 
that there is a microhystcrcsis loss (which is independent of the frequency) 
giving a decrement of 0,01 . The initial slope of the 5 t vs/ curve is then about 



